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We study the influence of entanglement on the relation between the statistical entropy of an open 
quantum system and the heat exchanged with a low temperature environment. A model of quantum 
Brownian motion of the Caldeira-Leggett type - for which a violation of the Clausius inequality has 
been stated by Th.M. Nieuwenhuizen and A.E. AUahverdyan [Phys. Rev. E 66, 036102 (2002)] - 
is reexamined and the results of the cited work are put into perspective. In order to address the 
problem from an information theoretical viewpoint a model of two coupled Brownian oscillators is 
formulated that can also be viewed as a continuum version of a two-qubit system. The influence of 
an additional internal coupling parameter on heat and entropy changes is described and the findings 
are compared to the case of a single Brownian particle. 



I. INTRODUCTION 

Open systems are subject to dissipation of energy and 
fluctuations in their degrees of freedom. Within the the- 
ory of quantum dissipative systems P, Q the starting 
point in describing noise and damping is the Hamiltonian 
fftot = Hs + He + HsE where the Hamiltonian of the to- 
tal system -fftot is expressed as a sum of the Hamiltonian 
of the subsystem of interest Hs, a Hamiltonian He mod- 
eling the environmental degrees of freedom and an inter- 
action term Hse- For quantum objects this coupling to 
the environment leads in addition to the phenomena of 
decoherence and entanglement. In the low temperature 
respectively strong coupling regime these entanglement 
effects become important. Under the unitary evolution 
of the density operator ptot(0 = U{t,0)ptot{0)U^ {t,0) 
with U{t,0) ~ ex-p[~ jrHtott] an initial product state of 
subsystem S and bath E evolves into a correlated state 
with ptot(i) 7^ Psit) ® PE{t) for t > 0. For zero temper- 
ature the closed, total system is in its ground state and 
therefore the von Neumann entropy S{ptotit)) stays zero 
for all times. But for t > the pure state of the whole 
system is an entangled state of subsystem and bath with 
S{ptot{t)) < S{psit)) + S{pE{t)). The subsystem there- 
fore is in a mixed state with S{ps{t)) > even for zero 
bath temperature. 

In most applications, especially in quantum optics 
0, 01 , the coupling between system and bath can be as- 
sumed to be weak which allows neglecting entanglement 
effects (Born-Markov- Approximation) and applying the 
formalism of Markovian quantum master equations. In 
this case the statistical thermodynamics of the open sub- 
system are governed by the quantum Gibbs distribution. 
In the strong-coupling quantum regime, the stationary 
state of the subsystem ps{t) is still Gaussian but non- 
Gibbsian due to the entanglement with the bath. If the 
total Hamiltonian is harmonic, a Gaussian initial state 
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of the subsystem remains Gaussian for all times. Its 
density matrix ps{t) is completely characterized by the 
first and second moments of the relevant operators. The 
Heisenberg equation of motion for these operators is the 
quantum Langevin equation [^. The characterizing 
moments are determined by the stationary values of the 
quantum Langevin equation and can be calculated alter- 
natively by applying the quantum fluctuation-dissipation 
theorem. The statistical entropy associated with that 
stationary quantum state is the von Neumann entropy. 
The systems exchange of heat with the environment ~ 
which is defined as the change in energy due to redistri- 
butions in phase space - is related to the thermodynamic 
entropy by the Clausius inequality. This thermodynamic 
entropy can only be identified with the statistical entropy 
when PS takes the form of the canonical density matrix. 
This is just the case for negligible interaction between 
subsystem and environment. 

From an information-theoretical point of view above 
considerations become important. The Landauer princi- 
ple 0, II] which is based on the Clausius inequality states 
that "many-to-one" operations like erasure of informa- 
tion require the dissipation of energy. Deleting one bit 
of information is accompanied by a released amount of 
heat of at least fcTln2. This erasure is connected with 
a reduction of entropy, and thus cannot be realized in a 
closed system. Therefore the information-carrying sys- 
tem has to be coupled to its environment. To avoid a 
rapid destruction of the necessary quantum coherence the 
quantum subsystem should be placed in a low tempera- 
ture environment. Thus the coupling might be relatively 
strong compared to thermal energy. Since the Landauer 
principle is dealing with information processing and era- 
sure, the relevant entropy is the statistical entropy of the 
system. Statistical entropy and heat are defined sepa- 
rately. Thus, the relation between both the quantities 
can be examined. 

The purpose of our paper is to study deviations from the 
Clausius inequality and Landauer bound respectively in 
the strong coupling quantum regime. An analytic treat- 
ment of this issue is given within the framework of the 
Caldeira-Leggett model of quantum Brownian motion. 
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In the first part of this paper we want to discuss 
quite controversial recent work ^| and put 

some of those results into perspective. Therefore the 
quantum Langcvin equation of a harmonically bound 
quantum particle which is based on the Caldeira-Leggett 
model is presented. The stationary moments which are 
obtained from this equation characterize the reduced 
density matrix which is used to define thermodynamic 
quantities. Then, changes in heat and in statistical 
entropy for adiabatic parameter variations are com- 
pared and the applicability of the Clausius inequality 
in the strong coupling/ low temperature quantum 
regime is discussed. In the second main part we will 
focus on a model of Brownian motion of two coupled 
oscillators, which can be understood as a continuum 
version of a two-qubit system [l 3j . With regard to 
recent work done on continuous variable computing 
|14L [Tsj l the impact on quantum information theory is 
studied. It will become clear that additional internal 
degrees of freedom can lead to different results as in 
the case of quantum motion of a single Brownian particle. 



II. CALDEIRA-LEGGETT MODEL OF 
QUANTUM BROWNIAN MOTION 

Brownian motion is a prominent example of an open 
quantum system (16| . The standard model of quan- 
tum Brownian motion, often refereed to as the Caldeira- 
Leggett model 0,01, is a system-plus- reservoir model. 
The whole system is governed by the Hamiltonian 
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where q and p are the Heisenberg operators for coor- 
dinate and momentum of the Brownian particle mov- 
ing in the harmonic potential V{q) — ^muif^q^ . The 
particle is coupled to a bath of N harmonic oscillators 
with variables Xi and pi and uniformly spaced modes 
LVi = zA. The interaction is bilinear in coordinates of 
system q and bath Xi. For the couplin g-pa rameters Ci 
the so-called Drude-UUersma spectrum [2 with large 
cutoff-frequency T and coupling constant 7 is chosen: 
Ci = ^2^miUjfl^T'^/i:{ijjf + F^). The bath is character- 
ized by its spectral density J{ij->), which takes the form of 
the Drude spectrum J(w) = ^ujT'^ / [uJ^ + F^) in the ther- 
modynamic limit (sending A — s- 0, iV — > 00 and keeping 
F ~ iVA =const.). The potential renormalization term 
J2 / {'^^i^i)q'^ ensures that V{q) remains the bare po- 
tential. Neglecting this self-interaction term, the positive 
definiteness of the total Hamiltonian Htot would just be 
guaranteed for 7 < mwp /F and - since F is large - would 
restrict the applicability of the model to weak-coupling 
approximations (7 <^ mujo). 



A. Quantum Langevin equation 

From the Hamiltonian the Heisenberg equations of 
motion for the operators q and p and the bath variables 
Xi, Pi are received. By eliminating the bath degrees of 
freedom the quantum Langevin equation 0| of a par- 
ticle moving in the potential V{q) can be derived: 

mq{t) + ^^+J^ dt'lit ~ t' W) - ry(t)-g(0)7(t). (2) 

The stochastic character of this integro-differential equa- 
tion with the friction kernel ^{t) = 7Fe~'"'*' comes into 
play by considering the initial distribution of the bath 
variables which determines the noise term rjit): 

w / fO) \ 
■q{t) = a I Xi{0) cos Wit H — sinujit I . (3) 
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Assuming an uncorrelated initial state with the reservoir 
being in canonical equilibrium at temperature T — /3~^, 
Pe exp(— ri(t) is a stationary Gaussian operator 
noise with {ri{t))pj^ = and the correlation function i3||: 



K{t ~ t') 
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In the case of initial correlations of particle and bath 
the correlation function contains additional terms which 
affect the dynamics on the timescale t < 1/F. To fully 
characterize the reduced dynamics it is thus important 
to specify the initial preparation. 



B. Stationary state 

The relaxation dynamics of the moments {q^{t)) and 
{p^it)) described by eq. end up in a stationary state 
for t ^ 00. The stationary correlations can be calcu- 
lated alternatively by applying the quantum fluctuation- 
dissipation theorem which establishes a connection 
between the quantum mechanical dynamical susceptibil- 
ity x{ijj) = x(t-t')e"^* = [mwg - mw^ - iujj{uj)] 
and the equilibrium fluctuations (g^) and (p^): 
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If the dissipative part of the susceptibility x"{^) of the 
non-Markovian damped oscillator with three characteris- 
tic frequencies Ai, A2 and A3 (poles in the complex plane) 
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is inserted, an analytic expression for eq. (0) and © is 
derived 0: 
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where ipix) is the Digamma function and Aq = A3, 
A4 = Ai. The stationary state of the Brownian parti- 
cle is fully characterized by the variances {Tj) and © 
which determine the stationary density matrix ps of the 
subsystem [Tll2]|: 



ps{q, q) = 
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This reduced density matrix is different from the canoni- 
cal equilibrium density matrix pth ^ exp(— /JiJs) for any 
finite coupling 7. The statistical entropy of the quantum 
- the von Neumann entropy S{ps) - is jol. |23|: 



state Ps 

S{ps) = -Tr[ps Inpg] = - y^p„ Inpn = 



(10) 



with Boltzmann constant set to fc^ = 1 and the subsys- 
tems phase space volume v defined by 
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as well as the eigenvalues of ps{q^ q'), 

= l/{v + 1/2) [{v - l/2)/(« -I- 1/2)]" , (12) 

which are obtained as solution of the problem 
/ dx' ps{q,q')fn{q') = Pnfniq), whcrc the eigenfunctions 
/„ are given by /„ ~ ^/cHn{cq)e~'^ ' ^^fiHi Hermite 
polynomials and c — [(jp^) / {h^ {q^))Y ^ ^ ■ 

The von Neumann entropy of the subsystem is in- 
creased with raising the coupling 7 at a given temper- 
ature (fig. Even at T ^ the subsystems entropy is 
larger than zero. This effect is due to the correlations be- 
tween subsystem and bath which prevent the subsystem 
from reaching a pure state for T — > 0. The probabilities 
to find the subsystem in an exited state depend on the 
coupling to the environment ,23'|. In the weak coupling 
limit, where ps = pth, expression H1U|) gives the entropy 
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of an harmonic oscillator in canonical equilibirum. 

It is important to remark here, that the entropy H1U|) 

deviates from the difference of the total entropy S'(ptot) 
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FIG. 1: Temperature dependence of the entropy expressions 
S{ps) a,nd Sp for different values of the system-bath cou- 
plings 7 (in units of muiQ/V). Dark lines from bottom to top: 
Sips)f^o; S{ps)-f=i; S{ps)-,=5- Grey lines: S'p,^=i (dashed) 
and Sp^-y^s. Other parameters: luq = 1, rn — 1, T — 10, h — 1. 



and the entropy of the bath in absence of the particle 



S{pE)'r=o which is given by [2j 
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In the same way the thermodynamic potentials Fp and 
Up can be derived which are related by Sp — l3{Up — Fp). 
One can see from fig. ^ that Sp vanishes at T — > 
whereas S{ps) does not. Proceeding in this way, the 
entropy Sp also contains the part of entropy that is asso- 
ciated with the quantum mechanical correlations or par- 
ticle and bath. Since this conditional entropy is nega- 
tive for entangled systems, the statistical entropy of the 
Brownian oscillator alone is underestimated by Sp. Thus, 
in our further treatment we will concentrate on the en- 
tropy S{ps)- 



C. Thermodynamics of adiabatic changes 

Nieuwenhuizen and Allahverdyan examined the 
validity of the Clausius inequality respectively the Lan- 
dauer bound in the strong coupling quantum regime. 
They found a violation of these principles at very low 
temperatures due to the existing correlations between 
subsystem and bath. With respect to quantum infor- 
mation theory they concluded that quantum mechani- 
cal information carrier therefore could be more efficient 
than their classical counterparts. A controversial subject 
in this context is the appropriate choice of heat and en- 
tropy expressions. Our purpose in this section is to clarify 
this issue and to put the findings of Nieuwenhuizen and 
Allahverdyan into perspective. 

The internal energy of the Brownian oscillator can be 
defined as the mean energy in the stationary state 



Us = Tr [Hsps] = (Hs) 
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This expression differs from the equivalent to H14() defined 
internal energy Up. The difference Up — Us can be inter- 
preted as the interaction energy Uint (7^ {Hse)^-) which 
is related to the free energy Fp by 
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Choosing the parameter values m, loq and T as in fig. 
n the ratio r — Ui-at/Up at zero bath temperature is 
given by r(^=,„^2/r) « 0.03 and T(^=5„<^2/r) ~ 0.10. 
For kT — Hluq the ratios are Tf^-y^muj^/r) ^ 0-01 and 
'r(7=5m^;5/r) ~ 0.05 respectively. 

The total differential dUs of the internal energy Us, 

dUs = TT[psdHs] + Tr[Hsdps] ^ SW + SQ (17) 

can be divided into two parts [25|. The first term re- 
sults from the change of the parameters m and luq in the 
Hamiltonian, so it is a mechanical, non-statistical object 
and will be referred to as work SW: 



5W = muJo{Q^)duJo 



(18) 



The second term TT[Hsdps] represents the variation of 
Us due to the statistical redistribution of the phase space, 
which will be associated with the change in heat 6Q: 

SQ = S^Q + S„iQ with 
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\2 am 2m dm J 
Now the validity of the Clausius inequality 
5Q < TdS 
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can be evaluated. The Second Law of thermodynamics 
in the formulation by Clausius states that in a quasi- 
static process, during which the system at all times 
passes through equilibrium states, one has dSth = SQ/T. 
The thermodynamic entropy Sth defined by the Clausius 
equality can only be identified with the statistical en- 
tropy S{ps) at thermal equilibrium where ps = Pth with 
pth = Z~-^ exp {—f3Hs) and Z — Trexp {—(3Hs), because 

dS = -Tr[(ipthlnpth] = Tr[dpthlnZ] + l3TT[dpthHs] = 
= /?Tr[pthHs] = f^SQ - dSth- (22) 

Concerning the Landauer principle, which is based on 
the Clausius inequality, but deals with information pro- 
cessing and erasure, the relevant entropy is the statistical 
entropy S{ps)- Therefore we want to compare changes 
in the statistical entropy dS to changes in heat SQ in- 
duced by adiabatic variation of the systems parameters. 
The total differential of the von Neumann entropy H10() 
is given by 



dS = — Tr[(ip5 Inps] — In 



dv. 



(23) 



Thus, the sign of the change in S{ps) is determined by 
the sign of the change in v. Note here that the param- 
eters TO and loq are chosen as independent quantities, so 
that we can examine SQm < TdSm and SQ^^ < TdS^j 
separately. 
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FIG. 2: Changes in heat SQm (positive) and entropy-term 
TdSm (negative) versus bath temperature T for different val- 
ues of the system-bath coupling 7 (in units of mcoQ/T): 0.5 
(dashed), 1 (grey), 5 (dark). The oscillator parameters are 
chosen to be uio — 1 and m — 1. The cutoff-frequency is set 
to F = 10 and ft = 1. 



Fig. [3 shows the temperature dependence of the changes 
in heat 6Qm and of the term TdSm for different coupling 
strength 7. While dQ„i is always positive, which means 
that the Brownian particle absorbs heat during an adia- 
batic increase of mass, the change of entropy and there- 
fore the product TdSm remains negative. In the high- 
temperature limit one has SQm — > as well as TdSm — > 
and therefore the behavior of an uncoupled harmonic os- 
cillator characterized by (|13() . The term TdSm converges 
relatively slowly towards zero because of the increasing 
factor T. The smaller the coupling 7, the faster is the 
convergence of the two terms SQm and TdSm ■ If the tem- 
perature T goes to zero, then the product TdSm does 
as well. The amount of heat SQm exchanged with the 
bath stays positive even in this limit and equals —dmUint- 
Thus, the Brownian particle can extract heat from the 
bath even at T = O^a fact that was already extensively 
discussed in ref . la [l3 ■ 



Fig. Ogives results of the cited work Nieuwen- 
huizen and AUahverdyan studied the influence of 
adiabatic changes in mass on the phase space volume v 
given by eq. (|ll|l . At T = one receives a monotonous 
decreasing function which should converge with in- 
creasing mass to the minimal value Vmin — 1/2 which 
results from the uncertainty relation. They found, that 
at moderate temperatures the phase space volume first 
decreases for low masses, then reaches a minimum and 
finally increases nearly linearly with the mass. The 
increasing phase space volume means a positive sign for 
the entropy change dS. Therefore the authors conclude 
that the different signs of SQ and TdS would only occur 
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FIG. 3: Phase space volume v versus mass m for different 
T-values. From bottom to top: T = 0; 0.1; 0.2; 0.25. Oscil- 
lator potential V{q) = |ag^ witli fixed spring constant a — 1. 
Other parameters: 7 = 1, F = 500, h = 1. (See also ref. ) 
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FIG. 4: Phase space volume v versus oscillator mass m for 
different temperature values. From bottom to top: T — 
0; 0.25; 0.5; 1. Oscillator potential V{q) = \mLolq^ with 
frequency uoq — 1. Other parameters: 7 = 1, F = 500, h — 1. 



at very low temperatures due to quantum correlations 
between system and bath. 

In contrast to that, fig. ^ shows the mass dependence 
of the phase space volume using the moments defined 
by eq. (T)) and ijHJi. One can clearly see that even for 
moderate temperatures the phase space volume does 
not increase, but reaches a temperature-dependent limit 
value. This value is given by the phase space volume of 
an uncoupled harmonic oscillator in canonical equilib- 
rium: vth = ^ coth/3?iwo/2. Both, sending m ^ 00 or 
coupling 7 — > 0, finally leads to the standard case of the 
quantum Gibbs distribution. 

The differences between our findings and the results 
in which become obvious in the figures |21 and 0] can 
be explained as follows: in the paper by Nieuwenhuizen 
and AUahverdyan the harmonic potential V{q) = 
imcjgg^ is expressed by ^aq'^ with spring constant a. 
Varying m and keeping a fixed leads to the results of fig. 
121 But since the limit m ~* 00 and 7 — > should lead 
to the same result of w = 1/2 at T=0, this choice of the 
potential is inconsistent. In fig. |21the phase space volume 
at T = does not reach the value v = 1/2 even for high 
masses (instead v « 0.6 for the given parameter values). 
In order to receive the correct expressions for {q'^), (p^) 
and U in the weak coupling or high temperature limit 
one has to set a = mujQ , which has also been done in the 
cited work in different contexts. 

Choosing the potential V{q) — ^mujQq'^, which is then 
affected by the variation of the mass, shows that the 
anomaly of different signs is an even stronger effect than 
found in [T]| . This indicates that not only quantum cor- 
relations at low temperatures might play a role but also 
classical correlations between the damped oscillator and 
its environment at moderate temperatures. 



III. QUANTUM BROWNIAN MOTION OF 
TWO COUPLED HARMONIC OSCILLATORS 

In order to study the infiuence of additional param- 
eters on the results stated above we introduce a model 
of quantum Brownian motion of two coupled oscillators 
which can be viewed as a continuum version of a two 
qubit system. In contrast to former work concerning the 
relaxation dynamics of two coupled oscillators 13, 26] we 
focus on the stationary state. The Hamiltonian Hs' of 
the open quantum system S' now reads 

Hs' ^Ha + Hb + Hab, (24) 

where Ha and Hb are the Hamiltonians of the two har- 
monic oscillators A and B with masses majmh and fre- 
quencies uja and uj(,. Hab describes the interaction be- 
tween them. Before deriving a quantum Langevin equa- 
tion it is necessary to discuss different couplings between 
the oscillators as well as between the system Hs' and the 
bath He and to choose an appropriate model. 

A. Coupling between the two oscillators 

In the framework of quantum optics the coupling be- 
tween oscillators is often chosen to Hab = —Dqaqb with 
coupling parameter D. In this case one problem is the 
constraint D < ^ma'mhOJa'-Ob as a condition for real eigen- 
frequencies of the system which restricts the range of al- 
lowed parameter variations. In our further treatment we 
will concentrate on the interaction Hamiltonian 

HAB^^Diqa-qbf. (25) 

This Hamiltonian is clearly inspired by its mechanical 
analogy - a restoring force proportional to the relative 
distance of the two oscillators - and leads to real eigenfre- 
quencies of the system Hs' for all values of the coupling 
parameter D. 
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FIG. 5: Oscillators separately coupled to the bath: eigenfre- 
quencies v of eq. I28II for a finite model with A'' + 2 = 10 os- 
cillators versus the system-bath coupling strength j (in units 
of mV). The critical value 7crit is given by eq. 12911 . Param- 
eters: Ua = 2, ujb — 5, m = ma,b = 1, D — 2 and A = 1, 
r = iVA = 8. 
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FIG. 6: Bath-coupling to the center of mass: eigenfrequencies 
V of eq. (I31II for a finite model with N + 2 = 10 oscillators 
versus the system-bath coupling strength 7 (in units of mV). 
Parameters: uja ~ 2, cui, = 5, m — ma,b = 1, D = 2 and 
A = l, r = iVA = 8. 



B. Coupling between system and bath 

In order to study the case of strong system-bath- 
coupling 7, the positive definiteness of the Hamiltonian 
has to be guaranteed in the range of relevant 7-values. 
Therefore we will discuss different couplings between sys- 
tem and bath in the following section. 



1. Oscillators separately coupled to the bath 

If each oscillator is coupled separately to the bath ac- 
cording to the coupling in Q then we receive the follow- 
ing interaction Hamiltonian Hs'e'- 
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The equation for the eigenvalues v of the total system 



reads: 



Ha + He 
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Fig. El shows the influence of the system-bath-coupling 
strength 7 on the eigenfrequencies of a finite system con- 
sisting of the two oscillators coupled to an environment of 
eight oscillators. The lowest eigenvalue vi decreases with 
increasing coupling strength and becomes imaginary at 
a critical value 



(28) 



7crit 



yJimaOJl + D){mbUjl + D) - D 



(29) 



which means exponentially increasing amplitudes and 
therefore instability of the whole system. In the ther- 
modynamic limit this critical value becomes very small, 
so that this model is only suitable in the weak coupling 
case. 



2. Bath-coupling to the center of mass 

A more adequate model is the bath-coupling attached 
to the center of mass R of the system iJg/ . The system- 
bath interaction could then be described by the Hamil- 
tonian: 
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(30) 



This coupling leads to the following eigenvalue equation 
of the total system -fftot: 
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In Fig. I^one can recognize that the lowest eigenvalue is 
only slightly reduced by increasing the coupling strength 
and remains real for all 7-values. 
Since this interaction term assures the positive definite- 
ness of the total Hamiltonian we will use this system- 
bath coupling in the further examination. Additionally 
this coupling allows us to transform the system easily to 
normal coordinates which simplifies the analysis in the 
case of identical oscillators. 
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C. Langevin equation of two coupled Brownian 
oscillators 



By transforming the Haniiltonian Hs onto coordinates 
for the center of mass R — 1/M{maqa + mtgh) and the 
relative coordinate x = {qa — qb) (with total mass M and 
reduced mass fi) and eliminating the bath variables, the 
following system of coupled equations for the Heisenberg 
operators x and R can be written down: 

X = -^lx{t)-{ujl~uol)R{t) 



R = -nlR{t)-^J^ ^{t~t')R{t)dt' 



with the frequencies 



ruatol 



(32) 
7(0^(0) " 
(33) 

(34) 
(35) 



and damping term "f{t) and noise term ri{t) as defined 
in section II. A. Solving eq. (|32|l as an inhomogeneous 
differential equation and inserting the solution into l|33|) 
gives a Langevin equation for the center of mass R with 
new damping term ^(t) and new noise term fjlt): 



MR + 



dV(R) 
dR 



dt'^{t - t')R{t') = f){t) - i?(0)7(t) 

(36) 



which describes the motion of R in the effective potential 
V{R) = l/2M^\R'^ with frequency 



(37) 



influenced by the damping 



1^ ^^^^ ^f'^\ osn,it~t') (38) 



and the stochastic force 
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i=l 



Xi(0) cos(ajit) + 3^^-^ siTL{ujit) 
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x{Q) cos(17^i) + sin(f7^i) 



(39) 



In the case of identical oscillators the equations (|32|l and 
(|33|l are decoupled. The relative coordinate performs a 
harmonic oscillation and for R the Langevin equation of a 
Brownian particle with mass M and oscillator frequency 
is received, 



MR+Ainj^R' 



dt'-f{t - t')R{t') 



77(t)-i?(0)7(i), 
(40) 

which is equivalent to quantum Langevin equation Q in 
the first part of this paper. 



D. Stationary state 

In order to calculate the stationary correlations for the 
general case we again apply the quantum fluctuation- 
dissipation-theorem. 

From eq. (|36|l one obtains the dynamical susceptibility 
x'lii^) — [Mflj^ — Muj'^ — iuj^{uj)]^^ and can express the 
variance of the center of mass by 



(41) 



-yV^oj 1 

/ duj—^r—:: — t; — - COt\l( — BfUJJ^ 



where 



2\1 



(42) 



Fig. [7|shows the temperature dependence of the variance 
(|41|l in comparison to the limiting cases: 



(i?'>th 

(i?')cl 



— ^cothi/3nf7fl for 7->0 (43) 
2A/f2/j 2 



kT 



for T > 



(44) 



-0.06 



-0.02 




. 1 



0.2 0.3 



. 4 



. 5 



FIG. 7: Temperature dependence of (-R^) 14111 for parameter 
values uja = 2, oji, = Suja, ma = mb — 1, D = Milfj, F = 10 
and 7 = 5Aff2|{/F, (dark line). For increasing temperature 
(R^) converges to the limit cases {R^)th (grey line) and {-R^)ci 
(dashed line) as given by eq. 1431 1 and 1441 . 



In the same way we can specify the variance of the 
center of mass momentum Pr : 



h r°° 1 
(Pb) = TT^'^^ / rf^Xfl(t^)c^'coth(-/3fiw) 



(45) 



27r 
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as well as the variance of the relative coordinate x 

} f°° 1 
(^') = TT / dwhcoth{-f3hL,)x'^{uj) = (46) 

- r H (^g-^fc)" l^^UJ coth{^ /3huj) 
and the corresponding momentum 

(47) 

The stationary correlation {xR) is obtained by transform- 
ing on normal coordinates y — C,x + R, z — x + -dR. 
Because of {yz) = in the stationary state one receives 



{xR) = 



ax^)+HR') 

1 + 



(48) 



where C Mn'^_ ~ {mbul + rUaUjl + M / iiD) 



MVl\_ - {maujl + rubuj^) 



(49) 



and U,± are the normal frequency of system H24|l with 
Hab given by if^ : 



(free) oscillation of the relative coordinate so that the 
values (a;^) and {p^.) are determined by the quantum 
Gibbs distribution of an u ncoupled o scillator with mass 
/i and frequency = a/w^ + D / ji. The motion of R 
is described by eq. and leads to the stationary 

variances given by lO and 0. 

In this case of identical oscillators the von Neumann 
entropy S{ps') of the system Hs' can be expressed as 
sum of the entropy of the center of mass coordinate Sr 
and the entropy of the relative coordinate Sx ■ 



4^ ( Hi! + _L 

TO6 rua 



Further correlations, such as (RPr), {xp^), {Rpx), ■■■ are 
zero in the stationary state. 



S{ps') = Sx + Sr, 



(53) 



where Sr and Sx are defined similar to equation Hl()|l 
with the phase space volumes 



'"x = \/{x'^){pl)/h and 



(54) 
(55) 



The exchange of heat 5Q and the change in entropy 
dS{ps') are defined equivalent to eq. lfT7|l and by 



SQ = Ti-[Hs'dps'] and 
dS = -Tr[(ips' Inps']- 



(56) 
(57) 



We now want to study deviations from the the Clausius 
inequality in the case of identical oscillators. Regarding 
variation of the mass M the Clausius inequality reads: 



(50) 



SQm < TdS_ 



M- 



(58) 



With regard to an information theoretical viewpoint, 
we use again the statistical entropy instead of the 
thermodynamic entropy S'th for which the Clausius 
equality - by definition - is fulfilled for quasi static 
processes. 



E. Thermodynamic of adiabatic changes 

The stationary Gaussian state of the subsystem is com- 
pletely characterized by the correlations 1)41(1 . (|45l) - (|48|l . 
The internal energy again is defined as the stationary 
mean value of the systems Hamiltonian 124|l : 

+ Uml{R^) + p{u;l - u;l){xR). (51) 

In case of identical oscillators (ma.b — M/2, uja.b = (^) 
the internal energy Us' turns into 

Us' = ^{P^} + lMcj^{R') + ^hn^coth{^f3hn^), (52) 

where {Pr) and (i?^) are given by the moments defined 
in equations iQ and © with oscillator parameters M 
and Lu. We can apply the weak coupling limit to the 



In the weak-coupling case of a single harmonic oscil- 
lator with mass M, entropy and heat are not affected 
by adiabatic variations of the mass: 6Qm = TdSn = 0. 
This is different in the case of two coupled oscillators 
weakly interacting with the bath, where the additional 
coupling parameter D leads to an increase of 5Qm 
and TdSM with rising temperature (narrow dashed 
curves in fig. andlHb). Nevertheless for 7 ^ the 
equality SQm = TdSM holds for all values of T and D. 
Furthermore, in this weak coupling approximation it is 
5Qm>Q and dSn > for aU T and D. 

The impact of a non-zero system-bath-coupling 7 
can be studied from fig. |SJl: the exchanged amount of 
heat 5Qm is increased at a given bath temperature T, 
whereas the term TdSM is reduced. At low temperatures 
the product TdSM becomes negative, so that changes in 
heat and entropy have different signs. This effect is the 
larger the stronger the system-bath-coupling 7 is. At 
high temperatures kT 3> ftw both terms become equal. 
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0.2 0.4 0.6 0.8 1 

kT/huo 



(a) Coupling parameter D between the oscillators 
A and B is set to D = Mu)'^. 



. 1 




0.2 0.4 0.6 0.8 1 

kT/huj 

(b) Weak coupling between the oscillators A and 
B with parameter D — O.OIMlj^. 

FIG. 8: Changes in heat SQm and entropy-term TdShi ver- 
sus bath temperature T for two different values of the oscil- 
lator coupling D. The system-bath coupling is chosen as in 
fig. m From top to bottom: SQm,-,^5, SQm,^=i, SQm,j^o.5, 
rdS'Af,7=o.5, rdS'Af,7=i, TdSM,'t=5- The oscillator parameters 
are oj = cua.t = 1 and M — 2ma.b ~ 1- The cutoff-frequency 
is set to r = 10 and h = 1. (Compare to fig. 



If the coupling D between the two osciUators is 
reduced (D 0), the system behaves hke a single 
Brownian particle as could be supposed by comparing 
the figures IHb and [21 

For D oo the system again behaves like a single 
Brownian oscillator. This is shown by fig. and fig. 

which give the D-dependence of 6Qm and TcLSm 
for different temperatures T and couplings 7. One 
can recognize that the values of 5Qm for D ^ 00 and 
D ^ Q are equal (as well as the values for TdSu) 
and correspond with the results for a single Brownian 
oscillator. 

Furthermore, fig. |^ shows that, depending on the 



.08 




^-0.0 



2 4 6 8 10 

D/Mlo'^ 

(a) 5Qm and TdSM versus coupling D at bath 
temperature T = hu/k. 



. 02 




0.5 1 1.5 2 2.5 

(b) SQm and TdSM versus coupling D at bath 
temperature T = 0.5huj/k. 

FIG. 9: Changes in heat 5Qm (grey line) and entropy-term 
TdSM (dark line) versus coupling parameter D at two differ- 
ent bath temperature T in comparison to the case of ^ 
where 5Qm,~,=o ~ TdSM, -y^o (dashed line). Parameters are 
chosen as in fig. |^ 



chosen parameter values of oscillators (M, cu) and bath 
(7, F, T), there may exist a range of ZJ-values where the 
terms TdSM and SQm have equal signs. Fig. ^jp makes 
clear that at lower temperatures such a region does not 
exist necessarily. 

Summarizing the results of this section, we can state 
the following: 

Compared to the case of a single Brownian particle, an 
open quantum system with internal degrees of freedom 
shows additional effects. In our introduced model the di- 
rection of heat and entropy fiow due to mass variations 
depends on the coupling strength between both the os- 
cillators. Already at moderate temperatures the flow of 
heat and entropy occurs in the same direction, whereas 
in the model of single Brownian motion this is reached 
only in the high temperature and weak coupling limit 
respectively. Therefore varying the coupling parameter 
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offers the possibility to adjust the ratio of heat exchange 
and change in the subsystems entropy. 
Of course, the resulting changes in heat and entropy de- 
pend on the chosen interaction between the oscillators. 
As pointed out at the beginning of the section, this cou- 
pling had to be selected carefully to ensure the positive 
definiteness of the total Hamiltonian. Nevertheless this 
model provides a good example for studying the influence 
of additional degrees of freedom. 

IV. SUMMARY AND CONCLUSIONS 

We have discussed the statistical thermodynamics of 
quantum Brownian motion of two coupled oscillators in 
the strong coupling quantum regime and have compared 
the results to the case of a single Brownian particle. 
In both cases quantum correlations between subsys- 
tem and bath lead to deviations from the canonical 
equilibrium thermodynamics. The quantum Langevin 
equation which has been derived for a system of two 
coupled Brownian oscillators describes the evolution 
of the Heisenberg operators. The stationary moments 
of these operators characterize the reduced density 



matrix completely. This density matrix contains all the 
accessible information about the quantum state. The 
statistical entropy of this state is measured by the von 
Neumann entropy. Reducing this entropy by quasi-static 
parameter variations is equivalent to a decrease in the 
information which is gained by a measurement. With 
regard to continuous variable quantum computing we 
examined the relation between changes in the subsystems 
statistical entropy and the exchange of heat with the 
environment. We found that this relation deviates from 
the Clausius (in)equality at low temperatures due to the 
existing correlations between system and bath. Related 
results of former work 0, ^| were nevertheless 
put into perspective. Concerning quantum information 
processing, the validity of the Landauer principle which 
is based on the Clausius inequality but deals with 
the statistical entropy seems indeed questionable - at 
least for open quantum systems which are non-weakly 
interaction with a low temperature environment. 
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